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Abstract 

We present a Monte Carlo (MC) study of and Zq production of lon- 
gitudinally polarized proton collisions at a center of mass (CMS) energy of 
500 GeV. All results were obtained with the SPHINX MC code. We consider 
two different types of polarized parton distributions with large and small gluon 
polarization. Various spin asymmetries are found to be sensitive to the shape 
of the polarized gluon distribution between 0.1 ~ x ~ 0.5. The asymmetries 
are approximately constant with respect to the transverse momentum carried 
by the and Zq, except if the first moment of the gluon polarization is 
large. In this case some asymmetries show a significant variation with 
transverse momentum. 
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I. INTRODUCTION 



In the past few years polarized deep inelastic scattering became a very attractive topic 
in high-energy physics. Experimentally, an important advance has been the measurement 
of the structure function g\ of the proton with high statistics [0,0. On the theoretical side 
much effort has been done to decompose g\ into parton distributions, but the situation is 
still unclear (see, e.g., || for a recent review). This is especially true for the polarization 
of the sea quarks, the polarization of the gluons, the shape of the parton distributions at 
intermediate Bj0rken x, and the small-x behavior of the parton distributions. 
There are in principle means to calculate all these quantities either from QCD sum rules or 
from lattice QCD. It will be exciting to compare future calculations and experiments. 
Meanwhile the parton distributions are constructed from models of the nucleon, e.g., from 
Carlitz-Kaur type models @] and many others. An even more heuristic way to obtain the 
parton distributions is simply to assume a generic shape of the polarized distributions and 
fit the strength and size parameters to the g\ data. This procedure is of course in no 
way unique. Many of these distributions are motivated from properties of the unpolarized 
distributions at larger x. Some fits invoke an SU(3) symmetric polarized sea in order to fit 
the first moments of the parton distributions in terms of /5-decay and hyperon decay matrix 
elements, see, e.g., ||. 

The naive way of comparing model parton distributions with parton distributions defined 
through expectation values of field theoretic operators has been much criticized ||,f?] . It has 
been argued that a proper comparison can only be made after the axial anomaly has been 
taken into account which shifts the first moments of the naive parton distributions against 
the first moments of the field theoretic parton distributions by 

A <?/ = A g/ - ^-Ag, (1) 

where Aqf denotes the field theoretic distribution of a quark with flavor f and Aqf denotes 
the corresponding quantity within a quark model. Ag is the first moment of the polarized 
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gluon distribution. Equation ([[]) is valid for any light quark flavor f. This equation is special 
for longitudinal polarization. Analyzing the data within a quark model, one has to fix the 
four unknown quantities Am, Ad, As, and Ag, neglecting the contribution of heavy quarks. 
But even when assuming an SU(3) symmetric sea quark distribution, the g\, /3-decay, and 
hyperon decay experiments only fix three out of the four unknown first moments. The 
shapes of the distributions are of course not known at all from these considerations. 
Therefore other experiments have been proposed to complement the measurement of g\ of 
the proton. The most important ones being semi-inclusive deep inelastic scattering || and 
polarized proton scattering • In the present work we concentrate on proton scattering. 
The production of large transverse momentum photons in polarized proton collisions might 
be a suitable observable, but there is a large background from pion decays in the final state, 
as was investigated in detail in ref. | IT| . This is not so for and Z production at RHIC 
where the background from other events will be very small ||. The limiting factor in this 
experiment will be the number of and Z that can be produced. There are at least 
three different poorly known properties of the partons, which must be entangled in future 
experiments, namely, the amount of flavor breaking in the sea quark distribution, the size 
of the gluon distribution, and the shape of all these distributions. In this paper we only 
consider two different parametrizations. The first parametrization stems from a Carlitz- 
Kaur type constituent quark model with a polarized gluon, radiatively generated sea quarks 
starting from Q 2 = 10 GeV 2 , and a slightly broken SU(2) x SU(3) symmetry 0. This model 
does not use the anomaly equation ([[]), but rather identifies Aqj with Aqj. This can be 
justified a posteriori because in this model the gluon polarization Ag = 0.262 turns out to 
be small. The interesting feature of these distributions is that they satisfy the Bj0rken sum 
rule and the Ellis- Jaffe sum rule. 

The other parton distributions have been taken from ref. ||. These distributions also have 
SU(3) symmetric sea quarks, but the sea quarks are somewhat stronger polarized due to 
the lower starting point Q 2 = 4 GeV 2 of the evolution. The analysis in ref. || proceeds 
via the anomaly equation (P and a larger ratio F/D = 0.590. Therefore, a much larger 



gluon polarization of Ag = 1.971 at Q 2 = 4GeV 2 is obtained. The basic question we want 
to discuss is: Can one distinguish between these parametrizations in a W ± /Z Q experiment, 
e.g., at RHIC? 

The measurement of the flavor breaking in the sea quark distribution by studying W ± and 
Z production is not addressed in this paper. A discussion of this interesting subject can be 
found in ref. fll^f . 

The paper is organized as follows. In section II we give a very brief review of the SPHINX 



MC code [pLl|,p^3| , in which the electroweak matrix elements have now been included. In 
section III we present the parton distributions that are used in our MC simulation. Finally, 
in section IV we present and discuss the results. 

II. THE MONTE CARLO PROGRAM SPHINX 

The application of perturbative QCD to Drell-Yan like reactions is complicated through 
the possible emission of soft and collinear radiation of partons in the initial and final state. 
The importance of these bremsstrahlung effects is seen most clearly from the experimental 
transverse momentum spectra of W ± and Zo production cross sections. These spectra cannot 
be explained through an intrinsic transverse momentum of the initial partons. We must 
rather think of the transverse momentum as being generated by bremsstrahlung effects. It 
is not possible to calculate the K-factor exactly to arbitrary order a s , so at some stage one 
must rely on approximations. 



The MC generator SPHINX [11], which is an extension of the well-known PYTHIA 5.6 MC 



code [14]], treats parton radiation in leading logarithmic approximation (LLA). This means 
that it resums the planar bremsstrahlung graphs in an axial gauge. These graphs are exactly 
the universal bremsstrahlung contributions that do not depend on the particular process, 
e.g., ud — > W + or uu — > 7/Z0. In SPHINX and PYTHIA the bremsstrahlung cascade is 
resummed in backwards direction, starting from the most virtual partons immediately before 
the reaction takes place. The longitudinal momenta of the mother partons are reconstructed 



4 



using first order polarized GLAP equations and evolved parton distributions. All information 
from the GLAP equation necessary to reconstruct the longitudinal momenta of the mother 
partons is contained in the inverse Sudakov-like formfactor 

giving the probability that the parton remains at x from i> = /^(Q^/A 2 ) to t < = ln(Q"^/A 2 ), 
where Q\ and Q> denote the virtualities of the daughter and mother partons respectively 
The indices a and h denote the flavors of the partons and h a and hi, their helicities. The 
polarized splitting functions P a ,h a ->b,h bl c,h c (z) can be found in ref. |15| . As it stands, the 
integrals diverge for the emission of soft radiation due to the poles in the splitting functions. 
To cure this, the soft region in the integral is cut out and treated separately. This causes 



no problem ||14| . Note, that SPHINX cannot evolve parton distributions but rather needs 
evolved partons as can be seen from equation (0). For the leading-order GLAP evolution 
we used the code developed by Gehrmann and Stirling || . 

The full four-momenta of the mother partons can be reconstructed recursively: Consider 
two partons with four-momenta Pi,P2> an d virtualities Q\ = —pi and Q\ = — pf, in their 
common CMS frame. Let Q\ > Q\. Then the momentum of the mother parton 3 of parton 
1 is reconstructed while parton 2 is not changed. Finally, parton 3 and parton 2 are boosted 
to their common CMS frame and the procedure starts again until both partons have a vir- 
tually smaller than 1 GeV. 

As is well known, the variables x used in the splitting functions and in the parton distribu- 
tions have no unique kinematical interpretation in LLA, see, e.g., fllfj. For technical reasons 



in SPHINX the same prescription as in PYTHIA is used, namely, the s-approach where 
s = X1X2S is the momentary CMS energy of to partons at any stage of the evolution and s 
denotes the overall CMS energy of the reaction. 

The last ingredients to be specified are the hard scattering cross sections. At low transverse 
momentum a weak boson is produced through a simple Breit-Wigner resonance. In SPHINX 
the treatment of the s-dependent width of the resonance is the same as in PYTHIA 5.6. We 



simply use that the polarized hard scattering cross sections are related to the unpolarized 
ones by 

a w ±(s, hi, h 2 ) = ~<W 2 -i <r%£*(s), (3) 

where hi, and h 2 are the helicities of the quark and antiquark respectively. A similar relation 
is true for Z production 



C 2 C 2 

cr Zo {s, hi, h 2 ) = S hlhat -! { S hl ,+i C 2 + L C 2 + ^^ d + Cl ^ ( 4 ) 



with Ci = g v + g a and Cr = g v — g a expressed through the axial and vector coupling 
constants g a , and g v respectively. If the transverse momenta of the vector bosons become 
large the collinear approximation in the treatment of initial state radiation becomes bad and 
one should use higher order matrix elements in combination with the radiation algorithm. 
For this purpose we implemented the next-to-leading order QCD and QED polarized cross 
sections in SPHINX. The QCD cross sections can be found in ref. |T7| but we also list them 
here: The annihilation process q(hi) + q(h 2 ) —> g + V ectorboson has to lowest order in ot s 
the scattering cross section 

at 9 s z \ t u tu ) 

The QCD Compton process q(h) +g(X) — > q + V ectorboson has in leading order a s the cross 
section 

da q{h)m _ 1 2 1 f (s , t 2M 2 u\ _,(M 2 -t 2 



i2 a s{ b-haY-^ + X h)[ - i + 1 + ^-2Xh^-^). (6) 

Here hi, h 2 , h(= ±1), and A(= ±1) denote the incoming parton helicities. For Z production 
one has to substitute 
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b u = - — ( —1 + -sin 2 6 w 



1 - ^sin*Bw} , (7) 



for a and b, where g = e/sinOy/- In the case of W ± production one has to substitute 



a u d = b ud = - Tr -^=cos9 c , 
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2y/2 



a us = b us = ~7^ sin ® c > ( 8 ) 
for a and b, where &c is the Cabbibo angle. 

III. THE PARTON DISTRIBUTIONS 

In this section we describe the main physics properties of the parton distributions that 
we used in the simulation. For details see refs. [|J and ||. These parton distributions have 
been slightly modified because we use different unpolarized parton distributions in both 
cases. We have taken those from Gliick et al. fl8fl . If the helicity distributions are required 



to be positive, the absolute value of the polarized distributions cannot be larger than the 
unpolarized ones, viz. 

f ± (x) = 1 -(f(x)±Af(x))>0 =>\Af(x)\<f(x) (f = q,g). (9) 

Whenever |A/(x)| gets larger than f(x) we substitute it by |A/(x)| = f(x). This modifies 
the parton distributions a little bit at large x as can be seen from fig. £| 
The first parametrization of the polarized parton distributions ||] is based on a Carlitz-Kaur 
type spin-dilution model for the proton's distribution functions. It is referred to as CKT. 
The important aspects of this model are: 

• The gluon distribution can be expressed through the valence quarks because of the 
spin-dilution idea. 

• It uses a non-standard ratio of the SU(3) coupling constants F/D = 0.49 ± 0.08. 

• The free parameters of the model are chosen to fulfill the Bj0rken and Ellis- Jaffe sum 
rule which turns out to be possible due to the above non-standard Fj D ratio. 
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• The fit to the gi(x, Q 2 ) data is made at Ql = 10 GeV 2 . There are no polarized strange 
quarks at this Q 2 . 

• The first moment of the gluon distribution Ag = 0.262 is rather small in this model. So 
it makes no difference if one uses the anomaly relation ([!]) or the naive identification. 

The shape of these parton distributions can best be read off from fig. [l]. They are shown 
in this figure for a low resolution of Q = 2 GeV as well as for a large resolution that is 
approximately equal to the mass. Note that we use rif = 3 for the number of flavors in 
the GLAP equation, whereas rif = 2 has been taken in ref. ||. 

The second type of parton distributions, taken from 0, is a simple fit to g\ using /3-decay 
and hyperon decay data to fix the first moments of the polarized distributions. For the 
shape of the distributions a generic function is assumed and the shape parameters are fitted 
to the gi data. The important aspects of this distributions are: 

• The sea is taken to be SU(3) flavor symmetric and vanishes at Qg = 4 GeV 2 . 

• The gluon polarization Ag = 1.971 is rather large, mainly due to the larger value 
F/D = 0.590 and the anomaly equation ([!]). 

• Three different shapes of the distributions are considered. They all fit the data equally 
well. These parton distributions are called SETA, SETB, and SETC. 

The shapes of the the distributions can be read off from fig. |2[ It can also be seen from 
this figure that the valence quarks in both parton distributions are quite similar in size and 
shape. 

IV. RESULTS 

Before showing our results, we have to fix some notation and introduce the spin asymme- 
tries. By (si, sz), si, S2 = ±, we will denote a spin configuration of the scattering protons. 
Si = + means that the proton which moves in the direction of the positive z-axis has a spin 
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that also points in the z-direction. si = — 1 means that the spin of this proton points in 
the negative z-direction. The spin of this particle is by convention equal to its helicity. The 
spin of the other proton is denoted by S2- Note that its helicity is always opposite to its 
spin, e.g., (+, +) means that the 'beam proton' has spin and helicity +, whereas the 'target 
proton' has spin + and helicity — . 

By a 2 — > 1 process we mean resonance production through a hard scattering cross sections 
and @. A 2 — > 2 process proceeds by definition through the annihilation amplitude 
(H), or Compton amplitude (|6]). In both caseswe generate higher order corrections through 
the MC algorithm explained in section [Hj. This means in particular that 2 — > 1 processes 
together with radiation contain those parts of 2 — > 2 processes that survive the LLA. There- 
fore, one should switch on either 2 — >• 1 processes or 2 — > 2 process, not both. For large 
transverse momenta 2 — > 2 processes must be used because they describe the emission of an 
additional parton more accurately. On the other hand, if the transverse momentum is low, 
not any event is accompanied by an additional parton radiation and one should switch to 
the 2 — » 1 processes. If we compare the cross sections, we find that the switch from 2 — > 1 
to 2 — > 2 processes should be made for transverse momenta of about pt = 4 GeV. The same 
Pt is obtained by comparing the asymmetries derived from 2 — > 1 and 2 — ^ 2 processes. 
After these preliminary remarks we define the asymmetries and show the results. The parity 
conserving asymmetry is defined through 

(+") + (-+) 



A PC 

LL 



We consider three different parity violating asymmetries, namely, 

(-+) 




(10) 



A PV2 
LL 



A PV3 
■"■LL 



(11) 



(12) 



(13) 



A test asymmetry 

A TEST = ( + + )-(--) (14) 
LL (+ + ) + (--) 

which vanishes trivially from rotation invariance, shows the quality of the calculation. 
The spin-asymmetries are plotted in figs. |3], and [5] against the transverse momentum of 
the vector bosons in the CMS frame of the reaction. For each orientation of the protons, 0.1 
million events have been produced []. This is a roughly realistic number for the RHIC spin 
experiment with a luminosity of 10 32 /cm 2 s and 1.5 years of running. Some asymmetries are 
quite sensitive to the structure functions. But one must be careful not to attribute this to the 
size of the polarized gluon distribution alone. First of all, one can see from the asymmetries 
corresponding to SETA, SETB, and SETC, which only differ by the shape of the gluon 
distribution, that the magnitude of the asymmetries is quite sensitive to the shape of Ag. See 
fig. ^ and note that all Gehrmann/Stirling distributions (full symbols) have the same Ag but 
give very different asymmetries. These asymmetries overlap with those obtained for the CKT 
distributions (open circles), which have much smaller Ag. Interestingly, there is another 
significant distinction between the distributions with small and large Ag: Some asymmetries 
derived from SETA, SETB, and SETC vary significantly with transverse momentum pt, 
whereas those derived from the CKT distributions stay approximately constant. One can 
see this most clearly from the parity violating asymmetries A^ 2 ' w and A^ 3,w . This 
effect can be traced back to the hard scattering cross sections eqns. ([|), and (|J). Even 
though the Compton process is suppressed by a factor of 3/8, it can become comparable 
to the annihilation process for large transverse momenta because the Compton scattering 
falls off less rapidly than the annihilation scattering. At very large transverse momenta 
the Compton induced asymmetry will saturate the spin asymmetry. Fig. [| shows that the 
transition region from an annihilation induced asymmetry to a Compton induced asymmetry 



This requires about 3 hours of CPU time for one spin orientation on an IBM-RISC6000 
workstation. 
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can be reached for the distribution SETB (and as well for SETA, and SETC). This is not 
possible for the CKT distribution, at least up to ~ 20 GeV. The transition can also be 
observed in the corresponding spin asymmetry fig. [7]. Finally, one can see from fig. [| that 
the range in x where the cross section is large lies approximately between 0.1 ~ x ~ 0.5 
with a strong maximum at about x = 0.1. 

V. CONCLUSIONS 

We have presented an extension of the MC generator SPHINX to polarized and Zq 
production. Using two different types of parton distributions in a MC simulation at y/s = 
500 GeV, we find that several details of the calculated asymmetries show an unexpected 
sensitivity to the polarized parton distributions. The variation of the spin asymmetries 
A^ 2,W and Ai£ 3,w with transverse momentum < 20 GeV are especially interesting. 
The upper limit in p? in our simulation is dictated by MC statistics. For p? > 4 GeV we 
use NLO hard scattering cross sections in combination with initial state radiation. The 
simulations show that not only the size but also the pt dependence of some asymmetries is 
sensitive to Ag. To make use of this information, one has to measure the total pt of the 
hadronic products in addition to the charged lepton from the W decay. 
We only showed results for the transverse momentum spectra because they are especially 
interesting for measuring the polarized gluon distribution. There many other important 
observables that can be generated by SPHINX. Rapidity distributions, for example, are 
sensitive to flavor asymmetries in the polarized sea quark distribution. 
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FIGURES 




FIG. 1. MC reconstruction of the CKT parton distribution with SPHINX at Q 2 = 4GeV 2 
(left), and Q 2 = 8100 GeV 2 (right). For comparison, also the unpolarized quantities are shown 
(solid lines). 
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FIG. 2. MC reconstruction of the parton distributions SETA, SETB, and SETC with SPHINX 
at Q 2 = 4GeV 2 (left), and Q 2 = 8100 GeV 2 (right). Solid lines show the SETA distributions, 
dashed lines SETB, and dash-dotted lines SETC. For comparison, also the unpolarized quantities 
are shown (solid line). 
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FIG. 8. The x distribution of VF^'s for unpolarized proton scattering at yfs = 500 GeV. Only 
events with pr > 5 GeV have been considered. 
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